Influence of Dzyaloshinskii-Moriya interactions on magnetic structure of a spin-1/2 

deformed kagome lattice antiferromagnet 
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Motivated by the recent neutron scattering experiment on Rb2Cu3SnFi2 [Nat. Phys. 6, 865 
(2010)], we investigate the effect of Dzyaloshinskii-Moriya interactions in a theoretical model for 
the magnetic structure of this material. Considering the valence bond solid ground state, which 
has a 12-site unit cell, we develop the bond operator mean-field theory. It is shown that the 
Dzyaloshinskii-Moriya interactions significantly modify the triplon dispersions around the V point 
and cause a shift of the spin gap (the minimum triplon gap) position from the K to V point in the 
first Brillouin zone. The spin gap is also evaluated in exact diagonalization studies on a 24-site 
cluster. We discuss a magnetic transition induced by the Dzyaloshinskii-Moriya interactions in the 
bond operator framework. Moreover, the magnetization process under external magnetic fields is 
studied within the exact diagonalization and strong coupling expansion approaches. We find that 
the results of all above approaches are consistent with the experimental findings. 
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INTRODUCTION 



Frustrated quantum magnets have been studied in- 
tensively to unveil possible novel phenomena which can 
be caused by the interplay of frustration and quantum 
effects. Spin-1/2 two dimensional (2D) antiferromag- 
net on the kagome lattice is a central example of such 
frustrated quantum magnets. On the theoretical front, 
various studies have suggested several candidate ground 
states for the magnetically disordered phase of the spin- 
1/2 antiferromagnetic Heisenberg model: a variety of 
spin liquids^— and valence bond solids^ - -®. Experi- 
mentally, a lot of effort has been made to realize the ma- 
terial system of such geometrically frustrated antiferro- 
magnet. The Herbertsmithitei^-SiiS&SI ZnCu 3 (OH) 6 Cl 2 
has the most ideal structure with uniform exchange cou- 
plings and shows no magnetic ordering down to much 
lower temperatures than the Curie- Weiss temperature. 
On the other hand, there may be several additional fac- 
tors that would affect the nature of the ground state. For 
example, the Dzyaloshinskii-Moriya interaction and im- 
purity spins may be present in the material^ - — Because 
of these complications, it has been difficult to pin down 
the true ground state of the system. These considera- 
tions also suggest that it is important to understand the 
effect of various perturbations to the ideal kagome lattice 
structure and magnetic anisotropies to identify the true 
ground state of the system. 

Rb2Cu3SnFi2 is another compound that was synthe- 
sized as an attempt to materialize the kagome lattice 
antiferromagnet. Unlike the Herbertsmithite, however, 
it has a deformed kagome lattice structure (see Fig. [TJ. 
Magnetic susceptibility and magnetization measurements 
show that it is a magnetically-disordered 2D spin-gapped 
system^ As a model to describe the system, a spin-1/2 
antifcrrmagnetic Heisenberg model was considered with 
four different exchange couplings due to its deformed 
structure, labeled by J±, J2, J3, and J4. The exchange 




FIG. 1. (Color online) The deformed kagome lattice with 
four different couplings for the Heisenberg and Dzyaloshinskii- 
Moriya interactions respectively. The DM vectors are as- 
sumed to be perpendicular to the lattice plane and their di- 
rection is along the positive direction of the z-axis when the 
orientation from i to j in Si x Sj is clockwise as indicated 
by the arrows in the lower right triangles. For the VBS state 
with the 12-site unit cell, the valence bonds are formed in 
the interaction links with the strongest couplings, J\ and Di 
(thick black). 



coupling constants are estimated by fitting the thermo- 
dynamic data: J\, </ 2 , and J3 are quite similar in mag- 
nitude (^200 K), while J4 is about half of the others 
(Ji > J2 > J3 3> J a)- The size of the spin gap is esti- 
mated to be A/fc B =21 K. 

The recent neutron scattering experiment in conjunc- 
tion with a dimer series expansion study has shown that 
the ground state of Rb2Cu3SnFi 2 is indeed the valence 
bond solid (VBS) with a 12-site unit cell^ The valence 
bonds arc formed in the strongest interaction links with 
the coupling constant J\ (black line segments in Fig. [IJ. 
Further, it turned out that the Dzyaloshinskii-Moriya 
(DM) interactions are not negligible in Rb2Cu3SnFi2 
as seen in the triplon dispersions. If there are only 
(isotropic) Heisenberg interactions, the system has SO (3) 
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global spin rotation symmetry so that its triplon dis- 
persions should be triply-degenerate. When there are 
(anisotropic) DM interactions, the degeneracy of triplon 
dispersions is lifted. It was observed in neutron scattering 
experiment that the triplon dispersions of Rb2Cu3SnFi2 
are only doubly-degenerate. This fact indicates that 
there exist significant DM interactions in Rb2Cu3SnFi2. 

The valence bond solid state with a 12-site unit cell was 
studied by some of us for the spin-1/2 antiferromagnetic 
Heisenberg model on the deformed kagome latticed In 
that study, the weakest interaction was controlled with 
J 4 = aJ and J\ — Ji = J3 = J. The VBS state with 
a 12-site unit cell turned out to have lower energy than 
the VBS state with a 36-site unit cell, which was sug- 
gested for uniform kagome lattice antiferromagnet (a = 1 
case), when a is less than 0.97 within the bond operator 
mean-field theory. Concerning Rb2Cu3SnFi2, the triplon 
dispersions were predicted as well. The predicted triplon 
energy bands are triply-degenerate due to the SO (3) sym- 
metry of the Heisenberg Hamiltonian and spin gap is 
found at the K point. However, the spin gap (the mini- 
mum triplon gap) was observed at T point in the recent 
neutron scattering experiment 4i This presents another 
evidence for the importance of the DM interaction in 
Rb2Cu3SnFi2. 

In this paper, we theoretically investigate the effect 
of the Dzyaloshinskii-Moriya interaction in the deformed 
kagome lattice antiferromanget via various methods. The 
model Hamiltonian for the system consists of the Heigen- 
berg and Dzyaloshinskii-Moriya (TLd) interactions 

between the nearest neighbors; 

Hj + n D = J2 J iA • S j ■ + J2 D ' S < x S r (!) 

(id) (id) 

For the Heisenberg interactions, we consider four differ- 
ent exchange couplings, J\ > J2 > J3 3> J a, depending 
on the interaction link as shown in Fig. [TJ As for the 
DM interactions, we use the result from series expan- 
sion study^ According to the series expansion, the DM 
vectors have small in-plane components compared 
to the out-of-plane component, and the overall behavior 
of excitation spectrum is determined by the out-of-plane 
component whereas the in-plane ones make minor correc- 
tions. Due to this reason, the DM vectors are assumed 
to have only the out-of-plane component in our study 
(Dij//z). Like the Heisenberg interactions, we consider 
four different DM interactions, the magnitude of which 
is labeled by D n (n = 1, • • • ,4). D n is assumed to be 
proportional to J„, i.e. D n /J n = d z (n = 1, • • • ,4). 
The strength of the DM interactions is controlled by the 
parameter d z in this paper. The direction of is deter- 
mined by the Moriya's rule^^r— to be along the positive 
z-direction if the orientation from the site i to j in x Sj 
is clockwise in each triangle of the kagome lattice. In the 
lower right part of Fig. [I] the orientations from i to j 
for the positive DM vector are denoted with arrows in 
two possible types of triangles in the kagome lattice. It 
must be noted that the model Hamiltonian has SO(2) 



global spin rotation symmetry in the xy-plane because 
every DM vector is along the z-direction. Therefore, the 
generator of spin rotation in the xy-plane, JV S^z, is a 
good quantum number of the Hamiltonian. 

To analyze this Hamiltonian, we take three different 
approaches: bond operator mean-field theory, exact di- 
agonalization, and strong coupling expansion. First, the 
bond operator mean-field theory is used to investigate 
the effect of the DM interaction on the triplon disper- 
sions. Within the bond operator theory, we consider the 
VBS state with the 12-site unit cell in Fig. [2] as the 
ground state and obtain the triplon dispersions. These 
disperions are shown to be consistent with the results of 
the neutron scattering experiment. The dispersions are 
twofold degenerate, which results from the SO (2) sym- 
metry of the Hamiltonian. By increasing the strength of 
the DM interactions (d z ) from zero, we observe that the 
position of the spin gap (the minimum triplon gap) in 
the first Brillouin zone shifts from the K to T point at 
d z = 0.09. Beyond d z = 0.09, the size of the spin gap 
decreases gradually, but the position of spin gap does not 
change from the K point. It is interesting to note that the 
spin gap is linearly decreasing as a function of d z after 
d z = 0.09. By performing linear fitting, we estimate the 
critical value c? ZjC = 0.381, at which the spin gap becomes 
zero. The value of d z in Rb2Cu3SnFi 2 was estimated to 
be 0.18 in the series expansion study. At d z = 0.18, the 
spin gap in the bond operator theory is 0.464 Ji. This 
value is about four times larger than the neutron scatter- 
ing result, 0.126 J\. But, it is not surprising because the 
bond operator mean-filed theory tends to overestimate 
the spin gap. 

Secondly, we conduct the exact diagonalization study 
for finite size clusters in order to calculate the spin gap 
more accurately. In the 24-site cluster, the spin gap is 
found to be 0.229 Ji, which is not far from the experimen- 
tal value. As a matter of fact, this value is the best among 
the results of various approaches in our study. Via the ex- 
act diagonalization, we also study magnetization process 
because magnetically-disordered spin-gapped systems of- 
ten show interesting behaviors in magnetization process. 
Exact diagonalization is restricted to small size clusters 
due to the large unit cell of the deformed kagome lattice. 
To overcome this limitation, we also employ strong cou- 
pling expansion approach. Strong coupling expansion is a 
degenerate perturbation theory which provides an effec- 
tive Hamiltonian near the magnetic fields strong enough 
to close the spin gap. This approach allows the larger 
system sizes than the clusters used in the exact diagonal- 
ization study because the full Hilbert space is projected 
into the subspace of the states relevant to the magnetiza- 
tion process. From the exact diagonalization and strong 
coupling expansion studies, we obtain qualitatively the 
same magnetization curves. The magnetization curves 
increase continuously from zero to saturated magnetiza- 
tion as the magnetic field is increased, which is consistent 
with the magnetization measurement. 

The outline of the rest of the paper is as follows. In 
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FIG. 2. (Color online) The dimer configuration for the VBS 
state with the 12-site unit cell. Dimers are denoted with the 
blue thick line segments. The 12-site unit cell contains six 
dimers labeled by 1, ■ • ■ ,6. In our convention, S_r in a dimer 
is located at the vertices of the hexagon denoted with circle. 



Sec. HIl we develop the bond operator mean-field theory 
for the deformed kagome lattice antiferromanget. Exact 
diagonalization study is presented in Sec. IIIII In Sec. IIV1 
we perform the strong coupling expansion. Lastly, results 
are compared with the recent experimental results and 
summarized in Sec. [V] 



II. BOND OPERATOR MEAN-FIELD THEORY 



In this section, we construct the bond operator mean- 
field theory for the VBS state with the 12-site unit cell. 
The dimer configuration for the VBS state is depicted 
in Fig. [2] where dimers are denoted with thick blue line 
segments and six dimers in the 12-site unit cell are labeled 
by ,6. 



A. Bond operator representation 

To introduce the basis for the bond operator represen- 
tation, we consider the following Hamiltonian of a pair of 
spins, Sl and S^, which corresponds to an isolated spin 
dimer with the coupling constants J\ and D\ in Fig. [2J 



'H-dimer — Jl&L ' SiJ 
1 



D x z • Si x Si 



1 



(Ji + iD^Sl.+Sr - + -(Ji - iD{)S L -SR t+ 
JiSl,qSr,o, (2) 



where S n ,± = S n>x ± iS rhV and 5„,o = S n , z (n = L,R). 
The energy eigenstates are given by 



S t|Q) 

tl\o) 

4 lo> 

tl |0) 



where 



|s) = -L( e W2 |n) - e- fa / 2 ||t)), (3a) 

= - in) , (3b) 

|to>-^(e ra/2 |n>+ e - W2 |lt», (3c) 
l*-> = IU>, (3d) 

(4) 



Ji + iD l 



The corresponding energies are 



etc - -\h 



1 



Jl+Dl 



(5a) 
(5b) 
(5c) 



It must be noted that |0) is not the pure spin sin- 
glet state and 4 |0) is not one of the pure spin triplet 
states unless a = (or Dx = 0). , t+, t^, and t_ 
are operators to create a bond particle in the associated 
state and satisfy the bosonic commutation rules. The 
subscript of t-operator, {+,0,—}, implies the quantum 
number of Sl,z + <S(r,z, which originates from the SO(2) 
spin rotation symmetry of the Hamiltonian in Eq. ([2]). In 
other words, {+,0,—} indicates the magnetic quantum 
number of the t-bosons. 

With the above bond operators, the spin operators, Sl 
and Sr, can be rewritten as follows; 



S l ,- 
Sl,q 
Sr,+ 
Sr,- 



1 

i 

71 



(s f i_ + t\s - t\t + 4*-) • e 
(sh+ + tls - V t+ + tlt ) ■ e 



-ia/2 



ia/2 



-^=(-sh- - t^s - t\.to + tit-) ■ e la ' 2 , 
V2 

4=(-* f *+ - t- s ~ 4*+ + t-to) ■ e- ia/2 , (6e) 
V2 



(6a) 
(6b) 
(6c) 
(6d) 



S 



R,0 



-(-sh - 4 s + t + t + - 



(6f) 



In the above representation, we set H = 1. There are 
two reasons why we choose the states ([3]) as the basis 
for the bond operator representation against the conven- 
tional one . 38 ' 39 First, if we use the conventional basis, 
the DM interactions with the coupling constant D\ do 
not appear in the quadratic part of the Hamiltonian. By 
taking the above basis, those DM interactions can be in- 
cluded from the beginning. Secondly, the ground state 
cannot be a pure spin singlet state due to the broken 
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spin rotation symmetry in the presence of the DM in- 
teractions. In the conventional representation, or when 
a = 0, i-bosons are sometimes called the triplons be- 
cause they represents the spin triplet states. Although 
all of the t-bosons in our representation are not in the 
pure spin triplet states, we will call them the triplons 
throughout the paper. 



B. Mean- field Hamiltonian 

In order to write the Hamiltonian ([1} in terms of the 
bond operators, it is helpful to arrange the Hamiltonian 
into the following form; 



(7) 



71=1 



where T-L n includes the interactions with coupling con- 
stants J n and D n . That is, 

^=EE [*S£(r)-SS(r) 

r m— 1 

+ £> 1 *.Sr(r)xSS(r)] ) (8a) 
%=EE [J 2 S£(r).S™ +1 (r) 

r m— 1 

+ £ 2 z.S£(r)xS™ +1 (r)], (8b) 
^ = E E [J 3 S£(r)-S£(r + R) 

r (m,n;R)e/ 3 

+ D 3 z .S£(r) x S2(r + R)] (8c) 
^ = EE [^S™ +1 W-Sr(r) 

r m— 1 

+ £ 4 5-S™ +1 (r)xS2 l (r)], (8d) 

where Sjj = S^,, and the set 1% is given by 

I 3 ={(1, 5; r c ), (3, 1; -r B ), (5, 3; r A ), 

(2,6;-r A ),(4,2;-r c ),(6,4;r B )}. (9) 

In Eq. © and ©, r implies the lattice vector, m and n 
are dimer indices within the unit cell, L and R are indices 
to distinguish spins in a given dimer, and 

ya = iax, 



r c =r B - r A , 



(10) 



where a is the lattice spacing. In our convention, spin 
R is located at the vertices of the hexagon denoted with 
circle in Fig. [5J Rewriting the Hamiltonian © in terms 



of the bond operators is straightforward. In addition, the 
following hardcore constraint is required at each dimer to 
restrict the Hilbert space into the physical one; 



sKs+ t\ti = l. 



(11) 



z=+,o, 



Moreover, the operators s and are replaced by a num- 
ber s to describe the VBS order, s denotes the conden- 
sate density of the s-bosons. Therefore, we consider the 
following Hamiltonian in the bond operator theory; 



H=Hj+Hd 

-EX> 



■* 2 + E Cwvw- 1 

l=+,0,- 



,(12) 



where r is the lattice vector, m dimer index within the 
unit cell, and fi the Lagrange multiplier for the hardcore 
constraint. In principle, the condensate density s and 
Lagrange multiplier \i can depend on r and m. How- 
ever, based on the symmetry of the underlying lattice, 
they are taken to be uniform on the lattice. Arranging 
the Hamiltonian with respect to the degree of i-boson 
operators, 



H = N uc e + 



quartic: 



where N„ r is the number of unit cells. 



e = 6 [e s s 2 + fi(l - s 2 )] 



(13) 



(14) 



and Hquad and H quar tic consist of quadratic and quar- 
tics terms of t-operators, respectively. Cubic terms are 
ignored in the above Hamiltonian because magnetically 
disordered VBS states are of interest here. 

The quadratic part of the Hamiltonian is considered 
first. After Fourier transformations, it can be written as 
follows. 

Hquad — HquadS + T~Lquad,2 + Hquad,3 + T~Lquad,i, (15) 

where 

6 

Hquad.l = E E E ( £t < ~ M)4m( k )^m(k), 
i= + ,0,- k m=l 

(16a) 

= E EE^+'^^wW 

;=+,0,- k m=l 

+ h.c. 

+ E EE^ +,, ^( k ) i U(- k ) 

i=+,0,- k m=l 

+ h.c, (16b) 
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1~tquad,3 

E E E Ws+i-iDs)^ 



/=+,0,- k (m,Ti;R)G/ 3 



•*L( k X„(k) 



E E E ta+i-u*) 



s V kR 



i=+,0,- k (m,ri;R)e/ 3 



/i.e., 



*L(k)*l,„(-k) 



(16c) 



E EE^+'-^" 2 



Z=+,0,- k m=l 



4,m+l( k )*i,m( k ) 



/l.C. 



E EE i( J 4+^ 4 ) 



Z=+,0,- k m=l 



*L+i(k)* t _ I . m (-k) 



/i.e. 



(16d) 



Every term in Eq. (fTB| conserves the z-component of the 
total spin, J^. Sj^, as the original Hamiltonian ([5]) does. 

The quartic part of the Hamiltonian is decoupled in 
such a way that the resulting mean-field Hamiltonian 
conserves Yli Readers interested in the details are 
referred to Appendix [X] The corresponding mean-field 
parameters are defined in (|A2[) . The resultant mean-field 
Hamiltonian can be arranged in the following form. 



Hmf = N uc {e + e PQ ) - J ]T ^TrM ; (k) 

i=+,0,- k 

+ ^ E E^M'W^k)' ( 17 ) 



i=+,0,- k 



where 



A,(k)= ki(k),--- ,« ; ,6(k),tl w (-k),-- - ,^ 6 (-k) 



(18) 



M/(k) is a 12x12 Hermitian matrix, and epQ is the 
mean-field contribution to the energy and defined in 
(|A4j) . At this point, it is important to note that the 
mean-field Hamiltonian describes the triplon excitations 
and it is block diagonalized into three triplon modes of 
/ = +,(),—. As mentioned earlier, / is the magnetic 
quantum number of the triplon, which reflects the SO(2) 
global spin rotation symmetry of the original and mean- 
field Hamiltonian. This fact can also be seen from the 



following expression of the magnetization; 



i 



E EE 



J.rf TO (k)t,, m (k). (19) 



k m— 1 



The mean-field Hamiltonian is diagonalized via the Bo- 
goliubov transformation- which transforms t-bosons to 
bosonic quasiparticles (7); 



r,(k) = T,(k)A t (k), 



(20) 



where 



r,(k)= [ 7 u(k),--- , 71,600,71,,! (-k),---, 7 ^ 6 (-k) 

(21) 

and T;(k) is the transformation matrix. The transfor- 
mation matrix satisfies the following relation due to the 
bosonic statistics of the quasiparticles. 



T,(k)I B Tj(k)=I fl) 



(22) 



where Is is the diagonal matrix with the first (second) six 
diagonal elements being +1 (—1). Because of this bosonic 
constraint, IsM/(k) is diagonalized instead of M/(k), so 
the eigenvalue problem has the following form; 



T z (k)I B M,(k)Ti- 1 (k) =Ifln,(k), 



(23) 



where 17; (k) is the diagonal matrix containing eigenval- 
ues wj jm (k) (m = 1, • • • ,6) in the first and second half 
of its diagonal. It must be noticed that Tj _1 (k) contains 
eigenvectors of the eigenvalue problem in its columns and 
the eigenvectors are normalized according to Eq. (|22|) . 
Diagonalizing the mean-field Hamiltonian with the above 
transformation leads to the following ground state energy 
per site. 



1 

12 



1 



4A„ 



J2 E TrM '( k ) 



i=+,0,- 



2N U 



E EE w, ^( k ) 



z=+,o,- 



k m=l 



(24) 



The energy is a function of the condensate density s, La- 
grange multiplier fj,, and mean- field parameters in (|A2I) . 
The ground state is determined by solving Eq. (|25[) and 
[2| self-consistently. 



de 



ds 



= 0, 



de a 



dfi 



= 0. 



(25) 



C. Ground state energy and triplon dispersion 

In this subsection, the results of the bond operator 
mean-field theory are presented. We use the coupling 
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-0.41 




FIG. 3. (Color online) Plots of (a) the ground state energy 
per site (e flr ) and (b) the spin gap (A) versus the strength 
of the DM interaction (d z ). In the plots, the color represents 
the position of the spin gap; light blue for the spin gap at the 
K point and red for the spin gap at the V point. The shift 
of the spin gap position from the K to F point happens at 
d z = 0.09. In figure (b), the dashed line indicates the linear 
fitting for the spin gaps occurring at the F point. According 
to the linear fitting, the spin gap is closed and a magnetic 
transition occurs at d z , c — 0.381. 



constants estimated by the series expansion study^I for 
Rb 2 Cu 3 SnFi2; 



J 2 = 0.95 Ji, 
J 3 = 0.85 Ji, 
Ji = 0.55 Ji, 
d z =D n /J n = 0.18, 



(26) 
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-0.448 


0.229 


2.809 


SC24 




0.47 


2.73 


SC48 




0.47 


2.73 


Neutron-series^ 




0.126 





where n ~ 1, • • • ,4. For the DM interactions, increasing 
d z from to 0.18, we investigate the effect of the DM 
interactions on the ground state energy and the triplon 
dispersions. 



TABLE I. The ground state energy per site (e S r), the spin 
gap (A), the critical magnetic field (h c ), and the saturation 
magnetic field (h 3 ) when d z — 0.18. 



The ground state energy is plotted for various values of 
d z in Fig. [3] (a) . As shown in the figure, the ground state 
energy decreases monotonically once we include the DM 
interactions. The energy at d z = is slightly off from 
the general tendency. This is due to the fact that we em- 
ployed different mean-field decouplings for the cases of 
d z = and d z ^ because each case has different sym- 
metry; SO (3) spin rotation symmetry when d z = and 
SO(2) symmetry when d z 7^ 0. How the mean-field de- 
coupling is done in each case is explained in the Appendix 
lAl From the plot of the ground state energy, it turns out 
that the DM interaction tends to stabilize the valence 
bond state with the 12-site unit cell in the deformed 
kagome lattice antiferromagnet. The ground state en- 
ergy per site at d z = 0.18 is —0.433 J\. 

Figure |4] shows the variation in the triplon dispersions 
due to the change in the strength of the DM interaction. 
When there is no DM interaction (d z = 0), the triplon 
dispersions with I = +,0,— are all degenerate due to 
SO(3) spin rotation symmetry. The spin gap happens at 
the K point in the first Brillouin zone. These behaviors 
are also consistent with the triplon dispersions with J\ = 
J2 = J3 = J and J4 = 0.5 J of Ref. 32. As we increase 
d z from zero, the spin rotation symmetry is reduced to 
SO(2) symmetry by the DM interaction, so the threefold 
degenerate dispersions are split into two parts (thick red 
and thin blue) as shown in Fig. |4] In fact, the triplon 
dispersion with / = (thin blue) is not affected by the 
existence of the DM interaction except for a small overall 
shift. For the others of Z = +, — (thick red), which are 
still degenerate, the degeneracies at the T and K points 
are lifted by the DM interactions. In particular, there is 
significant change around the T point so that the energy 
minimum (or the spin gap) shifts from the K to T point 
at d z — 0.09. The spin gap is plotted as a function of 
d z in Fig. [3] (b), where color denotes the position of 
the spin gap; light blue for the spin gap at the K point 
and red for the spin gap at the T point. As shown in 
the figure, the spin gap does not change much when it is 
located at the K point. However, after the spin gap shifts 
from the K to T point at d z — 0.09, it decreases linearly 
as a function of d z . The spin gap (A) at d z — 0.18 is 
0.464Ji. The triplon dispersions at d z = 0.18 in Fig. |4] 
show qualitatively the same behavior as the excitation 



d z =0 



d z =0.045 



d z =0.09 



d z =0.135 



d z =0.18 



a>(k)/Ji 
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FIG. 4. (Color online) Variation in the triplon dispersions caused by the DM interaction (d z ). Color indicates the magnetic 
quantum number of the triplon dispersion; thick red for I — +, — and thin blue for I = 0. When d z — 0, three modes with 
I = +, 0, — are all degenerate, so they are denoted with one color, red. The above plot shows that as the strength of the DM 
interaction increases, significant change occurs around the T point so that the energy minimum (or the spin gap) shifts from 
the K to T point at d z = 0.09. 



spectra obtained in the neutron scattering experiment; 31 
the spin gap at the T point and the lowest energy band 
being doubly degenerate. Thus, the triplon excitation 
spectrum is captured well in the bond operator mean- 
field theory. The above features can be obtained just 
from the quadratic Hamiltonian (|16[) without the quartic 
contributions. The quartic contributions generate only 
numerical corrections without changing the qualitative 
behaviors. 

Now we consider the cases with d z > 0.18. If we in- 
crease d z beyond 0.18, the energy gap becomes smaller 
and the lowest energy bands become linear around the 
r point as shown in Fig. [SJ Based on this fact, it is 
expected that the triplon dispersion becomes gapless at 
certain point of d ZtC . The critical value d z . c = 0.381 is 
estimated by fitting linearly the spin gap values beyond 
d z = 0.09 as shown in Fig. [3] (b). When the spin gap 
is closed, it leads to the condensation of the triplons at 



d z =0.18 d z =0.24 d 2 =0.33 
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FIG. 5. (Color online) Triplon dispersions beyond d z — 0.18. 
The lowest energy bands become linear around the T point as 
we increase d z . 



the K point and then causes a magnetic order in the sys- 
tem. It is the transition to a magnetically ordered state 
induced by the Dzyaloshinskii-Moriya interaction. 

Lastly, we comment on the behaviors of the VBS state 
and the triplon modes under the external magnetic field, 
h = hi. The VBS state has no magnetization, thus 
the change due to the magnetic field occurs only to the 
triplon modes. Under the magnetic field, the triplon 
modes with I = +, — are coupled to the magnetic field 
while the other one with I = is not affected as shown 
in the following expression of the Zeeman interaction; 

H h = -hz-Y^ (27) 

i 

= E EE(- r ^tJ k )^( k )- (28) 

Z= + ,0,- k m=l 

Now the triplon dispersions are modified as follows; 

(k) -> w Jjm (k) - I ■ h, (29) 

where w/, m (k) (I = +, 0, — ; m = 1, • • • ,6) is the triplon 
dispersion when h = 0. As we increase h, the lowest band 
with I = + mode hits the zero energy at the critical 
magnetic field h c = A, where A is the spin gap when 
h = 0. Beyond h c , the VBS state becomes unstable due 
to the gapless triplon mode with I = +, the condensation 
of which causes a magnetic ordering in the system. 

The critical magnetic field, or the spin gap, is 0.464Ji 
when d z = 0.18. Comparing with the experimental re- 
sult A exp = 0.126Ji)^i the critical field is overestimated 
in the bond operator mean-field theory. For more ac- 
curate estimation of the spin gap, we conduct the exact 
diagonalization and strong coupling expansion studies in 
the next sections. In those studies, we also investigate 
the magnetization process. 
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FIG. 6. (Color online) Clusters with 12 and 24 sites used for 
the exact diagonalization study. 



III. EXACT DIAGONALIZATION 

The deformed kagome lattice has the large unit cell of 
12 sites. Because of the large-size unit cell, it is difficult 
to study large cluster sizes for the exact diagonalization 
(ED) study. Thus, we are limited to small size clusters 
shown in Fig. [5] with periodic boundary condition. In 
each cluster, we diagonalize the Hamiltonian 



Hj + Hd+Hh, 



(30) 



each term of which was already introduced in Eq. (J7J, 
©, and ((23). Fig. [7]shows the full spectrum of the 12- 
site cluster when there is no magnetic field. Its ground 
state occurs in S z = sector with the energy per site of 
-0.453Ji, where S z = Yl^i^i,* an d N s is the number 
of sites. The spin gap, which is the energy difference 
between the ground state and the lowest energy state in 
S z = 1, is 0.371 Ji. These values are modified in 24-site 
cluster; the ground state energy per is -0.448Ji and the 
spin gap is 0.229Ji. It is interesting to note that the 
value of spin gap is not far from that estimated in the 
experiment, 0.126 J\, even though it is obtained in the 
small cluster of 24 sites. 

If we include the Zeeman interactions and change the 
magnetic field (h), the ground state selected by the sys- 
tem changes. The ground state is determined by the com- 
petition among the lowest energy states in the sectors of 
S z = 0, • • • , N s /2; we call these states \E (S Z )). When h 
is small, the lowest energy state of S z = 0, \E a (S z = 0)), 
remains as the ground state. The lowest one of S z = 
N s /2, \E (S Z — N s /2)), becomes the ground state when 
h becomes large enough. How the system moves through 
the lowest energy manifold {\E (S Z )) \S Z = 0, • • • , N s /2} 
in intermediate region of h determines the profile of mag- 
netization curve. The energies of those states, E Q (S Z ), 
consist of concave curve in E a vs. S z plot as guided by 
the dashed line in Fig. [7] In that case, the system goes 
through all of the lowest energy states 4^ As a result, the 
magnetization (M z = 2S Z /N S ) increases gradually from 
zero to one (see Fig. [8] (a)). This can be understood by 
noting that the ground state is determined by the follow- 
ing equation; 



If E (S Z ) is concave like that of 12-site cluster, then 
the minimum of E (S Z ) — h ■ S z moves gradually 
from S z = to S z = N s /2. If the energy 
curve of {\E (S Z )) \S Z — 0, - ■ ■ ,N s /2} were convex, the 
ground state would jump abruptly from \E Q {S Z = 0)) to 
\E (S Z = Ns/2)). The gradual increase of magnetization 
from zero is consistent with the experimental result^ 

Now, we discuss about the effect of the DM interaction 
on the magnetization profile. Fig. |8] (b) shows magneti- 
zations for two cases of d z = and 0.18 in 24-site cluster. 
A few changes are observed. First, the spin gap or the 
critical magnetic field (A = h c ) is significantly decreased 
by the DM interaction; from 0.350 J x to 0.229 J x . This 
fact is consistent with the result of the bond operator 
mean-field theory DM interaction increases the satura- 
tion magnetic field h s , at which the system becomes com- 
pletely magnetized, as shown in Fig. (5](b). As a result, 
it extends the range, h s — h c , where the magnetization 
process goes on. 



IV. STRONG COUPLING EXPANSION 

In this section, we derive an effective Hamiltonian 
when the magnetic field is strong enough to cause mag- 
netization in the system. It is used to investigate the 
magnetization process of the deformed kagome lattice 
antiferromagnet. The effective Hamiltonian is obtained 
by applying the degenerate perturbation theory to the 
original Hamiltonian. This approach is called the strong 
coupling expansion^ 

We first arrange the Hamiltonian (j5(I)) as follows; 



Hj+H D +H h =H + H', 



(32) 



1 12 66 220 495 792 924 792 495 220 66 12 1 



E/J, o 



\ 5 ! 

\ 



■ = / 

I I ' 



_d_ 

dSz 



[E (S Z ) -h-S z ]=Q 



(31) 
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FIG. 7. (Color online) Full energy spectrum of the 12-site 
cluster when there is no magnetic field. Dashed line provides 
a guide for the shape of the the function E (S Z )- 
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FIG. 9. (Color online) Lattice structure for the effective hard- 
core boson model in the strong coupling expansion study. The 
finite clusters enclosed by the red line are used to solve the 
hardcore boson model exactly. In the lattice, each point de- 
notes a dimer in the original lattice of spins in Fig. [2] The 
left cluster in the figure is called 24-site cluster and the right 
48-site cluster, referring to the number of sites in the original 
lattice. 
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The effective Hamiltonian is generated by applying the 
first-order perturbation theory on the degenerate sub- 
space V Q . The Hamiltonian can be mapped into the hard- 
core boson model. To this end, we employ 6-boson at 
each dimer and map the states \s) and \t+) into zero- 
boson and one-boson states as follows; 



FIG. 8. (Color online) Magnetization curves for the clusters 
with 12 and 24 sites. 



where 

Ho = Hj 1 + Hdi ~ Adimer *Si,2> (33) 

i 

4 

n' = J2 (Kjn +H Dn ) + (A dimer ~h)J2 Si,z- (34) 

n— 1 i 

In the above equations, Hj n is the Heisenberg interac- 
tion part with the coupling constant J n and Hd„ is the 
DM interaction part with D n . Adimer = e* + — £s , where 
e t+ and e s are from Eq. ([5]). H is regarded as the un- 
perturbed Hamiltonian, which provides proper basis to 
describe the magnetization process, and W is the per- 
turbation for H - Notice that H corresponds to the 
Hamiltonian of independent dimcrs under the magnetic 
field of h = Adimerj the eigenstates of which are given 
by direct product of dimer states in Eq. ©. The mag- 
netic field of h = Adimer in T~i is intended to make the 
states with e t+ and e s to have the degenerate ground 
state energy in each dimer. Then, H has the degen- 
erate ground subspace; V = | \tp) = |s),|i+)}. 
The other states are gapped out from the ground state 



&t=| t+ ) (s\, 
b=\s) (t + \, 
n=\t+)(t+\ 



(35a) 
(35b) 
(35c) 



where and b are operators to create and annihilate a 
boson respectively and n is the number operator for the 
bosons. For the hardcore bosons, the original lattice in 
Fig. [2] is transformed to the lattice shown in Fig. the 
dimers in the former are mapped into the points in the 
latter, so that the 12-site unit cell in the former becomes 
the 6-point hexagon in the latter. The first order effective 
Hamiltonian is written as follows in terms of the hardcore 
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bosons. 

= N uc ■ 6e s 

1 6 

r rn—1 

+ h.c. 

+ l(J 3+i D 3 )J2 E Hn(r)b n (r + H) 

r (m,n;R)G/ 3 

+ h.c. 

r m— 1 

+ /l.C. 

1 6 

+ 7J2 22 22 n m(^)n m +i(r) 

r m— 1 

+ -£ J s^2 X! f>m(fK(r + R) 

r (m,n:R)e/ 3 
1 6 

+ -J4^ E n. m+ i(r)n m (r) 

r m— 1 

6 

+ (A dimer -h)^2^2 n ™(r), (36) 

r m— 1 

where bj = b\ and n-; — n\. The hardcore boson model 
Hamiltonian consists of the hopping and nearest neigh- 
bor repulsion. In this effective Hamiltonian, the hopping 
amplitude is comparable to the repulsion energy as shown 
in Eq. (|36|). This fact makes it difficult to predict the 
ground state of the Hamiltonian. For this reason, the 
Hamiltonian is diagonalized exactly for two clusters as 
shown in Fig. [9j The cluster shown in the left of the fig- 
ure is the 24-site cluster and the 48-site cluster is shown 
in the right. These clusters are used for the exact diago- 
nalization studies of the effective Hamiltonian. Since we 
are studying the hardcore boson model, the total number 
of bosons is allowed to have some number between zero 
and the number of sites in the underlying cluster. To ob- 
tain the magnetization profile, we find the ground state 
at each magnetic field. In the hardcore boson model, 
the magnetic field controls the number of bosons like the 
chemical potential. If we increase the magnetic field, then 
the number of bosons in the ground state increases. In 
fact, the number of boson is nothing but the magneti- 
zation of the ground state as the following relation indi- 
cates; 

M « = 6j^EEM')- (37) 

r m— 1 

The magnetization of the 48-site cluster is plotted in 
Fig. [TU] (a) for several values of d z . Interestingly, the 
magnetization plateau at M z = 1/2 appears when the 
DM interaction is weak (d z — ~ 0.09). The magneti- 
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FIG. 10. (Color online) Magnetization curves obtained in 
the strong coupling expansion study. In figure (a), the mag- 
netization curve for each value of d z is shifted upward to avoid 
overlapping of the curves. Each curve starts from M z — at 
h = 0. 



zation plateau at M z — 1/2 becomes narrower and nar- 
rower and eventually vanishes as the strength of the DM 
interaction increases from d z = to 0.18. The same be- 
havior in magnetization is observed in the 24-site cluster 
as well. 

Although M z = 1/2 plateau is an interesting point of 
the strong coupling expansion study, we must be careful 
about the interpretation. First of all, the exact diago- 
nalization study does not show any magnetization step 
which can be considered as a plateau. On the other hand, 
the effective Hamiltonian in the strong coupling expan- 
sion is obtained within the first order perturbation the- 
ory, so there are many higher order interactions which 
are ignored. In this regard, higher order interactions are 
to be included in the effective Hamiltonian in order to 
describe the magnetization process in terms of the hard- 
core bosons more accurately. Nonetheless, the first order 
effective Hamiltonian captures the overall increasing be- 
havior of the magnetization curve, already found in the 
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exact diagonalization study (see Fig. [10] (b)). The ef- 
fect of the DM interaction on the magnetization is also 
consistent with the exact diagonalization result; the crit- 
ical field value h c is decreased and the saturation field 
strength h s is increased when we increase the strength 
of DM interaction. It turns out that h c = 0.47Ji and 
h s = 2.73Ji for both of the 24- and 48-site clusters at 
d z = 0.18. The critical and saturation field strengths 
from various calculations are listed in Table Q] 



V. DISCUSSION 

We investigated the effect of the DM interaction in 
a model appropriate for Rb 2 Cu3SnF 12 via various theo- 
retical approaches. Assuming the DM vectors are per- 
pendicular to the lattice plane, we consider the valence 
bond solid state with the 12-site unit cell as the ground 
state of the system, which was confirmed by the neutron 
scattering experiment. First, using the bond operator 
mean-field theory, we study the effect of the DM interac- 
tions on the VBS state with the 12-site unit cell and the 
excitation spectrum. It is shown that the ground state 
becomes more stabilized by the DM interactions. The 
DM interactions have significant effect on the excitation 
spectrum or the triplon dispersions. First of all, the DM 
interaction breaks the SO (3) spin symmetry down to the 
SO(2) spin rotation within the lattice plane. As a result, 
the triplon modes with / = +, — , which are still degener- 
ate, are separated from the I = mode, with significant 
modification of the spectrum around the T point by the 
DM interaction. This change induces the shift of the spin 
gap (the minimum triplon gap) position from the K to 
r point in the first Brillouin zone at d z = 0.09. The 
spin gap at the T point and shapes of the lowest three 
triplon dispersions at d z — 0.18 in Fig. [4] are character- 
istic features observed in the neutron scattering experi- 
ment with Rb2Cu3SnFi2- Therefore, the bond operator 
theory correctly captures the DM interaction effect in 
Rb 2 Cu 3 SnF 12 . 

In the deformed kagome lattice antiferromagnet, we 
can think about two kinds of transitions induced by the 
Dzyaloshinkii-Moriya interaction and magnetic field re- 
spectively. For the former transition, the transition point 
is predicted to be d z _ c = 0.381 by fitting linearly the 
spin gaps for d z bigger than 0.09 in the bond opera- 
tor theory. At the transition point, the gapless triplon 
mode with I = +, — are condensed to cause the magnetic 
order. At this transition, the SO(2) spin symmetry is 
spontaneously broken. In fact, there were several studies 
on such kind of transition and possible magnetic order- 
ing for spin-1/2 antiferromagnet on the uniform kagome 
lattice^ - — Among those studies, the exact diagonaliza- 
tion study showed that the critical value for the mag- 
nitude of the DM vector is D/J = 0.1, where J (D) 
is the coupling constant of the Heisenberg (DM) interac- 
tions in the uniform kagome lattice. 35 Beyond the critical 
value, the system enters the Neel ordered state. Compar- 



ing our results with their study, the magnetic transition 
induced by the Dzyaloshinkii-Moriya interaction in the 
deformed kagome lattice antiferromagnet seems to have 
similar structure but with a larger value for the criti- 
cal DM interaction. However, it should be noted that 
the bond operator theory tends to overestimate the spin 
gap. Comparing with the neutron scattering experiment 
result, the spin gap value from the bond operator theory 
is about four times larger than the experimental result 
(see Table [J). Considering the linearly decreasing behav- 
ior of the spin gap with respect to the DM interaction in 
the deformed kagome lattice antiferromagnet, the critical 
value d ZjC is also expected to be overestimated. The crit- 
ical point can be estimated to be d z . c = 0.23 by using the 
following fact and assumption; (i) the fact that the bond 
operator theory prediction for the spin gap is four times 
larger than the experimental value at d z = 0.18, and 
(ii) the assumption that the spin gap of Rb 2 Cu3SnFi 2 
is a linear function of d z and it has the same slope 
that is predicted in the bond operator theory. Accord- 
ing to this estimation, Rb 2 Cu3SnFi 2 is located pretty 
close to the critical point. This situation is very simi- 
lar to one of the Herbertsmithite ZnCu 3 (OH)6Cl 2 . The 
strength of DM interaction in the Herbertsmithite was 
measured to be D/J^ 0.08 via electron spin resonance 
experiment fSZ which is pretty close to its critical value 
D/J = 0.1 from exact diagonalization study^ In this re- 
gard, Rb 2 Cu3SnFi 2 and ZnCu3(OH) 6 Cl 2 have common 
feature of being close to the quantum phase transition 
induced by the Dzyaloshinskii-Moriya interactions even 
though they have different lattice symmetries. 

The other transition occurs via the magnetization pro- 
cess. Under magnetic field perpendicular to the lattice 
plane, the lowest I = + triplon mode becomes gapless 
at the critical magnetic field h c = A. Then, it is con- 
densed to induce the magnetization along the direction 
perpendicular to the lattice plane. The entire magnetiza- 
tion process is studied by both exact diagonalization and 
strong coupling expansion. In particular, the exact diag- 
onalization is used to estimate more accurate magnitude 
of the spin gap. 

From the exact diagonaliztion study, the spin gap was 
computed to be 0.229 Ji for the 24-site cluster with the 
periodic boundary condition. This value is not far from 
the experimental result, 0.126Ji. The system exhibits 
continuously increasing magnetization curve without any 
magnetization plateau regardless of the DM interaction 
strength. The DM interaction reduces the critical field 
value h c and increases the saturation field strength h s . 
As a result, it extends the range, h s — h c , where magne- 
tization process goes on. 

In the strong coupling expansion, we restrict the 
Hilbert space to the subspace of relavent states for 
the magnetization process and then derive an effective 
Hamiltonian by applying the first order perturbation the- 
ory to the subspace. The effective Hamiltonian is mapped 
to the hardcore boson model, which consists of the hop- 
pings and repulsions of the hadrcore bosons under "chem- 



12 



ical potential" given by the magnetic field. Exact calcula- 
tions on several finite clusters at d z = 0.18 show qualita- 
tively the same magnetization curve compared to the ex- 
act diagonalization result of the full Hamiltonian. Inter- 
estingly, M z = 1/2 magnetization plateau was observed. 
It is concluded, however, that the effective Hamiltonian 
has to include higher order interactions to describe the 
magnetization process accurately. 
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Appendix A: Mean-field decoupling 

In this appendix, we explain the mean-field decoupling 
of the quartic part of the Hamiltonian (fT2"|) . The quartic 
part can be arranged in following form. 



H quartic 
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lowing mean-field parameters; 
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TABLE II. (m,n), R„ 



ED- 



where (m, n), R m n, 4 m ii, and -B m „ are listed in Table [Til 

We decouple the quartic terms in such a way that 
the resultant mean-field Hamiltonian preserves the SO(2) 
spin rotation symmetry. To this end, we define the fol- 



The above parameters are defined for (m, n) listed in 
Table [TT] and assumed to be independent of the lattice 
vector r. With these mean-field parameters, the quartic 
part of the Hamiltonian is decoupled into the following 



13 



form. 



and 
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A few comments on the mean-field decoupling are in 
order. First, it must be noted that there is another pos- 
sibility in decoupling the last terms with B mn in (|A1I) . 
It is to decouple them by using the direct channel; 



*t , _+t , 



ft i _+t / 



— > 



(tt ,t+,i) - (*!,«-,<) 



i 1 ,t_ 



(A4) 



(A6) 

where i and j mean (r, m) and (r + R mn , n), respectively. 
However, the direct channel mean- fields do not make any 
effect because they cancel each other so that their con- 
tribution becomes zero. The cancellation originates from 
the time reversal symmetry of the original Hamiltonian 
([7]). Due to the symmetry, (t\ i t+ i i) = (t_ and the 

direct channel mean- fields have no effect. Secondly, the 
above mean-field decoupling is intended for the Hamilto- 
nian with nonzero DM interactions (d z ^ 0). For the case 
without the DM interactions (d z = 0), reader is referred 
to Ref. 32. 
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